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ABSTRACT 

We discuss the antiGZK effect in the diffusive propagation of ultra high energy 
protons in intergalactic magnetic fields, which consists in a jump-like increase of 
the maximum distance from which ultra high energy protons can reach an ob- 
server. The position of this jump, Ej « 2 x 10 18 eV, is determined exclusively 
by energy losses (transition from adiabatic to pair-production energy losses) and 
it is independent of the diffusion parameters. The diffuse spectrum presents a 
low-energy steepening approximately at this energy, which is very close to the 
position of the second knee observed in the cosmic ray spectrum. The dip, seen 
in the universal spectrum cLS cL SI enature of the interaction with the cosmic mi- 
crowave background radiation, is also present in the case of diffusive propagation 
in magnetic fields. 

Subject headings: UHE Cosmic rays, diffusive propagation, GZK cutoff. 

1. Introduction 

The GZK cutoff (Greisen (1966), Zatsepin and Kuzmin (1966)) is a steepening of the 
ultra high energy (UHE) protons spectrum due to the interaction with the cosmic microwave 
background (CMB) radiation. The presence of an intergalactic magnetic field can modify 
the GZK cutoff up to its absence in the case of very strong magnetic fields, Sigl et al 
(2004), Yoshiguchi et al (2003), (for a physical explanation of this effect see Aloisio and 
Berezinsky (2004)). The proton propagation in magnetic field can affect the observed UHE 
proton spectrum also at energies (much) lower than the GZK cutoff. The crucial parameter 
which determines the modification of the spectrum is the distance d between sources. If this 
distance is much less than all propagation distances, such as energy-attenuation length, Z att , 
and diffusion length Z difT , the spectrum is not distorted and has a universal (standard) shape 
(Aloisio and Berezinsky (2004)). This statement has the status of a theorem. 
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All these effects depend strongly on the strength of the large-scale intergalactic magnetic 
field (IMF), the knowledge of which still remains poor. The modes of the UHE-proton prop- 
agation vary between quasi-rectilinear propagation in a weak field and diffusive propagation 
in a strong magnetic field. The experimental data on IMF and the models of origin of these 
fields do not allow at present to choose even between the two extreme propagation regimes 
mentioned above. 

The most reliable observations of the intergalactic magnetic field are based on the Fara- 
day rotation of the polarized radio emission (for reviews see Kronberg (1994), Valle (1997), 
Carilli and Taylor (2002)). The upper limit on the Faraday rotation measure (RM) in 
the extragalactic magnetic field, obtained from the observations of distant quasars, gives 
RM < 5 rad/m 2 . It implies an upper limit on the extragalactic magnetic field on each as- 
sumed scale of coherence length (Kronberg (1994), Valle (1997), Ryu et al. (1998)). For 
example, according to Blasi et al. (1999a), for an inhomogeneous universe B lc < 4 nG on a 
scale of coherence l c = 50 Mpc. 

According to observations of the Faraday rotations the extragalactic magnetic field is 
strongest, or order of 1 /iG, in clusters of galaxies and radiolobes of radiogalaxies (Valle 

(1997) , Kronberg (1994), Carilli and Taylor (2002)). The largest scale in both structures 
reaches l c ~ 1 Mpc. Most probably various structures of the universe differ dramatically by 
magnetic fields, with very weak field in voids and much stronger in the filaments (Ryu et al. 

(1998) ). Superclusters seem to be too young for the regular magnetic field to be formed in 
these structures on a large scale l c ~ 10 Mpc. 

In the case of a hierarchical magnetic field structures in the universe, UHE protons with 
E > 4 x 10 19 eV can propagate in a quasi-rectilinear regime. Scattering of UHE protons 
occurs mostly in galaxy clusters, radiolobes and filaments. Deflections of UHE protons can be 
large for some directions and small for the others. The universe looks like a leaky, worm-holed 
box, and correlation with the sources can be observable (see Tinyakov and Tkachev (2001), 
where correlations of ultra high energy cosmic rays (UHECR) with BLLacs are found). Such 
a picture has been suggested by Berezinsky et al. (2002b). 

A promising theoretical tool to predict the IMF in large scale structures is given by 
magneto-hydrodynamic (MHD) simulations. The main uncertainty in these simulations is 
related to the assumptions concerning the seed magnetic field. 

The MHD simulations of Sigl et al. (2004) and Sigl et al. (2003) favor a hierarchical 
structure with strong magnetic fields. Assuming an inhomogeneous seed magnetic field 
generated by cosmic shocks through the Biermann battery mechanism, the authors obtain 
~ 100 nG magnetic field in filaments and ~ 1 nG in voids. In some cases they consider IMF 
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up to a few micro Gauss as allowed. In these simulations UHECR are characterized by large 
deflection angles, of the order of 20°, at energies up to E ~ 10 20 eV (Sigl et al. (2003), Sigl 
et al. (2004)). Thus, the scenario that emerges in these simulations seems to exclude the 
possibility of an UHECR astronomy. These simulations have some ambiguity related to the 
choice of magnetic field at the position of the observer (Sigl et al. (2003), Sigl et al. (2004)). 
The authors consider two strong local magnetic field B ~ 100 nG and a weak field 

B <C 100 nG. The different assumptions about the local magnetic field strongly affects the 
conclusions about UHECR spectrum and anisotropy. 

The essential step forward in MHD simulations has been made by Dolag et al. (2003). 
In this work the Local Universe is simulated with the observed density and velocity field. 
This eliminates the ambiguity for the local magnetic field, that is found to be weak. The seed 
magnetic field, used in this simulation, is normalized by the observed magnetic field in rich 
clusters of galaxies. The results of these constrained simulations indicate a weak magnetic 
fields in the universe of the order of 0.1 nG in typical filaments and of 0.01 nG in voids. The 
strong large-scale magnetic field, B ~ 10 3 nG, exists in clusters of galaxies, which, however, 
occupy insignificant volume of the universe. The picture that emerges from the simulations 
of Dolag at el. (2003) favors a hierarchical magnetic field structure characterized by weak 
magnetic fields. UHE protons with E > 4 x 10 19 eV can propagate in a quasi-rectilinear 
regime, with the expected deflection angles being very small < 1°. 

The case of strong magnetic fields up to 1 fiG has been studied in Sigl et al. (1999), 
Lemoine et al. (1999), Stanev (2000), Harari et al. (2002), Yoshiguchi et al. (2003), Deligny 
et al. (2003). The interesting features found in these calculations are small-scale clustering 
of UHE particles as observed by Hayashida et al. (1996), Hayashida et al. (1999), Uchiori 
et al. (2000), Glushkov and Pravdin (2001), and absence of the GZK cutoff in the diffusive 
propagation, when the magnetic field is very strong. Many aspects of the diffusion of UHECR 
have been studied in numerical simulation by Casse et al. (2002). 

The small-scale clustering allows to estimate the space density of the sources (Dubovsky 
et al (2000) and Fodor and Katz (2000)). The recent Monte Carlo simulations (Yoshiguchi 
et al (2003), Blasi and De Marco (2004) and Kachelriefi and Semikoz (2004)) favor a number 
density of the sources n s ~ (1 — 3) x 10~ 5 Mpc~ 3 with rather large uncertainties (Blasi and 
De Marco (2004)). 

Diffusive propagation of extragalactic UHECR has been studied already in earlier work. 
The stationary diffusion from Virgo cluster was considered by Wdowczyk and Wolfendale 
(1979), Giller et al. (1980) and non-stationary diffusion from a nearby source was studied by 
Berezinsky et al. (1990a), Blasi and Olinto (1999b) using the Syrovatsky solution (Syrovatskii 
(1959)) of the diffusion equation. In this case the GZK cutoff may be absent. 
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A very interesting phenomenon, caused by propagation of UHE protons in the extra- 
galactic magnetic fields, has been recently found by Lemoine (2004). It consists in a low- 
energy steepening of the spectrum of UHE protons at energies below 1 x 10 18 eV produced by 
a large diffusive propagation time (exceeding the age of the universe) to the nearby sources. 
In this paper, we shall discuss the anti-GZK effect in diffusive propagation of UHE protons 
which is responsible for this low-energy steepening and discuss the transition from galactic 
to extragalactic cosmic rays. In our calculations we shall follow, like Lemoine (2004), the 
theoretical approach of Aloisio and Berezinsky (2004). 



2. Diffusive propagation in the analytic approach 

The analysis below is based on the analytical solution of the diffusion equation, found 
by Syrovatskii (1959). Using a distribution of sources on a lattice, the diffuse flux can be 
calculated as the sum over the fluxes from the discrete sources i: 

where b(E) = dE/dt is the proton energy loss, summation goes over all lattice vertexes, L p 
is the proton luminosity of a source, q{E g ) = E~ la is the generation function, and K{^ g ) is 
the normalization coefficient equal to 7 9 — 2 if 7 9 > 2 and 1/ ln(_E max /_E min ), if 7 9 = 2 (all 
energies are measured in GeV), and 

f Eg Die) 

is the Syrovatsky variable, which has the physical meaning of the squared distance traversed 
by a proton in the observer direction, while its energy diminishes from E g to E. From Eq. (1) 
one can see that the sources at distances r > 2^J\{E, E g ) give negligible contribution to the 
flux. 

In our calculations we shall use also the second Syrovatsky variable, which can be 
understood as the time needed by a proton to diminish its energy from E g to E: 

The Syrovatsky solution formally includes all propagation times up to t — > oo and the 
generation energies are restricted from above only by the maximum acceleration energy E^ x 
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Fig. 1. — Proton energy losses. The continuous line represents the sum of pair production 
and photopion production energy losses, while the dashed line gives adiabatic energy losses 
(l/E)(dE/dt) = H . The label Ep shows the energy where pair-production and adiabatic 
energy losses are equal. 

that a source can provide. In our case the propagation time from a source at fixed distance 
r must be smaller than the age of the universe t , and due to this condition one more upper 
limit on the maximum generation energy E™ x emerges. This limit is given by the condition 
t(E,E 9 ) < t and results in E g Iiax (E) < E 9 (E,t ), which can be calculated also by evolving 
energy backward in time from E at t = to E g at t = t . 

The upper limit _E™ ax in Eq. (1) is then the minimum between the two quantities: 
E g (E,to) and the maximal acceleration energy -E^ax; 



At small energies E < 2 x 10 18 eV E 9 (E,t ) < E^ x , while at larger energies E g (E,t ) > 
E^ x . In the calculations below we will assume E^ x = 1 x 10 22 eV. 

The crucial quantity in the following discussion, the proton energy loss [3(E) = (1/ E)dE/dt, 
is shown in Fig. 1. Note the characteristic energy Ep pa 2 x 10 18 eV, where the pair-production 
energy losses (3 e + e -(E) reach the adiabatic energy losses. 

Using these energy losses we can calculate E g ria,x (E). The results are presented in Fig. 



E g nax (E) 



mm[E g (E,t ), E^ x ] . 



(4) 
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Fig. 2. — Maximum generation energy i?™ ax defined as min[E g (E, to), E^ x ], where £™ is 
the maximal acceleration energy and t is the age of the universe (see text). 

2. At low energies E g (E, to) increases due to adiabatic energy losses. At the end of this stage 
the increase becomes more sharp because at large time t the pair-production energy-losses set 
in. Finally at E ~ Ep E g (E,t ) abruptly increases up to practically by a jump. The 
jump factor is roughly given by exp(t /f), where f is the energy- loss time which diminishes 
as the energy rises with the backward time. This behavior of E g na,x (E) is responsible for the 
antiGZK effect, which will be discussed in the next Section. 

We shall specify now the diffusion coefficient D(E), which determines X(E,E g ) in Eq. 
(2). In the following discussion we shall also use the diffusion length definition as: ld(E) = 



We assume diffusion in a random magnetic field with a strength B on the maximum 
coherent length l c , denoting this magnetic configuration by (B , l c )- This assumption deter- 
mines the diffusion coefficient D(E) at the highest energies when the proton Larmor radius, 



3D(E)/c. 



r L (E) > l c : 



D(E) 



1 crj(E) 
3 l c 



(5) 
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At "low" energies, when r L (E) < l c we shall consider three cases: 
(i) The Kolmogorov diffusion coefficient 



ZMB) = k(^r, (6) 



3 V 'c / 

(ii) The Bohm diffusion coefficient 

D B (E) = \cr L {E), (7) 

(iii) An arbitrary case D(E) oc E a , with a = 2 for the extreme energy regime. 

In all cases we normalize the diffusion coefficient by (l/3)c/ c at = l c . The characteris- 
tic energy E c of the transition between the high energy and low energy regimes is determined 
by the condition tl(E) = l c and is 

B , ^ 0.93 x 10 - (^) eV. (8) 

The smooth transition between the low-energy and high-energy diffusion regimes is 
provided with the help of an interpolation formula for the diffusion length: 

l m (E) = A d + d^l (9) 

with A d = r L (E) for the Bohm diffusion and A d = Icij'Ljlc) 1 ^ for the Kolmogorov regime. 

For completeness we shall give also the numerical expression for the Larmor radius: 

E 1 nG 

r L (E) = 1.08 x 10 2 Mpc. (10) 

v ' 1 x 10 20 eV B 1 K ' 

At distances r < ldm(E), the fluxes from individual sources i are calculated in the rectilinear 
approximation, and the diffuse flux is given by 

xrect/px L P K (lg) \- q g en(E g (E,rJ) dE g (E,r) 
J ? ( E) = ^^\ rj JE— 

where dE g /dE is given in Berezinsky et al. (2002a). 



3. Anti-GZK cutoff 

In this Section we shall demonstrate that, in contrast to the GZK cutoff, increasing of the 
proton energy losses at energy E > 1 x 10 18 eV results, in the case of diffusive propagation, 
in an increase of the maximal distance from which protons can arrive. 
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Fig. 3. — Maximal distance r max (E) to the contributing sources as function of the observed 
energy E. Three merging curves in the left-low corner give r^ x (£) and the dash-dotted curve 
gives r^^E), which numerically is very close to the energy-attenuation length l at t(E) = 
[{l/cE)dE/dt)\~ 1 . Continuous line is for the Kolmogorov diffusion D(E) oc E 1/3 at E < E c , 
dashed line - the Bohm diffusion D(E) oc E, and dotted line - the extreme high-energy 
diffusion D(E) oc E 2 . In the right panel D(E) = const case is also shown. Two different 
configurations of the magnetic field are considered: B = 1 nG, l c — 1 Mpc (left panel) and 
B = 100 nG, l c = 1 Mpc (right panel). 



We shall calculate below \(E,E™ X ), which according to Eq. (2) gives r^ ax /4, where 
r max{E) is the maximal distance from which protons with the observed energy E can arrive, 
as it follows from Eq. (1): 

\(E,E™) = J e de-±±. (12) 

In two extreme limits, at low energies and high energies, \(E, £™ ax ) can be calculated 
analytically. 

Let us start from the low-energy case E Ep, when only adiabatic energy loss operates. 
Using D(E) oc E a we obtain from Eq. (12) 

D(E) 



\(E,E t 



max\ 



aHn 



E 



(13) 



^max f ounc [ f rom conc lition r(E,E" iax ) = t is E™ x = Eexp(H t ), which results in 



< i a f x (^) = 2('^) 1/2 (e^-l) 1/2 , (14) 



where according to the WMAP data (Spergel et al (2003)) H t w I. 
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In the extreme high-energy regime E > 3 x 10 20 eV = E(dE/dt) 1 p=* 4.1 x 10 7 yr 
does not depend on energy and from Eq. (12) we have 



max v 



E: 



ace 



E r 



d,ff (S) = ^ L (15) 



Consider now the intermediate energies, when E approaches 1 x 10 18 eV, but E g (E,t ) 
remains less then E n m 4 x 10 19 eV, where photopion production starts. One obtains in the 
case D(E) oc E a 



rr x (^)oc V / ^(^f M ) a , (16) 

where r ee ~ Pe+ e -- I n this case r^ x (E) grows fast with E due to the fast growth of E g (E, to) 
(see Fig. 2). 

When E approaches Ep 2 x 10 18 eV, the value of r m ^ x is determined by the energy 
interval between E c and where -D(S) oc -E 2 . _E™ a:E there grows by a jump to E^ x , and 

r m : f x also grows by a jump to the high energy asymptotic value given by Eq. (15). 

The accurate numerical calculations are displayed in Fig. 3 for two different magnetic 
field configurations (1 nG, 1 Mpc) and (100 nG, 1 Mpc), respectively. 

In a diffusive regime of propagation there is an additional upper limit for a distance to 
a source, which we shall refer to as the rectilinear maximal distance r v ^ x (E). It is defined 

as 

= { * (B,£r) lltT < 17 > 

[ ci it t > to- 

At small E r r ^ x (E) = cr(E, i?™ ax ) = (c/H ) \n{Ef^/E) is larger than ct and r^(E), 
as one can see from Fig. 3. At large E r^(E) is smaller than r^ x (E), and thus the 
rectilinear upper limit becomes restrictive. 

The Syrovatsky solution (1) does not include automatically the restriction due to r™^ x (E), 
because propagation time there varies from to oo. The restriction (17) must be imposed 
in Eq. (1) additionally. This restriction is valid also in the case without magnetic field and 
numerically it is very close to the attenuation length l att (E) = E(dE/dl)~ 1 , which describes 
the ordinary GZK cutoff. 

Fig. 3 illustrates the antiGZK effect which we discuss here. While the energy-attenuation 
length l att (E) = E(dE j 'dl)~ 1 (or maximal rectilinear distance rj^.) diminishes with energy 
E and has the sharp GZK steepening at E ~ 5 x 10 19 eV, the diffusive maximum distance 
r m&x(E) increases with energy and has a sharp jump at energy Ej ^ 2 x 10 18 eV. As we 
discussed above, this energy is determined entirely by energy losses and it does not depend 
on the diffusion parameters. 
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The growth of r^ x (E) depends on the diffusive regime, as it directly follows from 
Eq. (14). 

4. Results and discussion 

The maximum distance r max (E) determines the number of sources which in principle 
can contribute to the observed diffuse flux J P (E): the flux from the sources at distances r 
larger than r max is suppressed as exp(— r 2 /r^ ax ). But inside the sphere with radius r max the 
fluxes from the sources are suppressed by X(E,E g ), which is less than X(E, £'™ ax ) and by 
E g la (E). By this reason, the jump in r max does not produce a jump in the flux at energy 
Ej. The situation is different at E < 2 x 10 18 eV, where r max (E) suppresses the diffuse flux, 
restricting the number of contributing sources. 

In Figs 4,5,6 we present the calculated diffuse spectra using Eqs. (1) and (11), in the 
case of two configurations (B , l c ) and for different distances d between sources. 

In our calculations the sources are located in the vertexes of a lattice, and summation is 
performed within the volume limited by r max (E) as described in Section 3. In fact, only the 
rectilinear limit is introduced by hand, while rjjj^(£) at lower energies appears automatically. 

As was expected, the energy of the low-energy steepening E s is nearly the same for 
all magnetic configurations and approximately coincides with the cross-over of adiabatic 
and pair-production energy losses Ep, and with the position of jump Ej. In accordance 
with r max (E) given by Eq. (14), the flux below the low-energy cutoff is the largest for the 
Kolmogorov diffusion (or D=const regime) and the lowest for D(E) oc E 2 diffusion, with the 
Bohm diffusion between them. 

In the calculations for a reasonable magnetic field configuration with Bq = 1 nG and 
l c — 1 Mpc, we have used a separation between sources d = 30 Mpc and d = 50 Mpc, which 
corresponds to a source space density 3.7 x 10" 5 Mpc" 3 and 8.0 x 10" 6 Mpc -3 , respectively. 
As was discussed in the Introduction, the small- angle clustering favors a density n s ~ (1 — 
3) x 10~ 5 Mpc -3 with some uncertainties. In the case of strong magnetic field Bq = 100 nG 
we have used a larger separation d = 100 Mpc to improve the agreement with observations. 

In Figs 4 and 5 we show the spectra in the case B — 1 nG and l c — \ Mpc. The critical 
energy where the diffusion changes its regime is E c ~ 1 x 10 18 eV, and the diffusion length 
at this energy is /diff ~ 1 Mpc. The best fit to the observations is obtained for 7 9 = 2.7 . 
The energy of the steepening in both cases is E s ~ 1 x 10 18 eV. The source luminosities L p , 
needed to provide the observed flux are very high, if one assumes a power-law generation 
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Fig. 4. — Energy spectrum in the case of B — 1 nG, l c = 1 Mpc and for the diffusion 
regimes: Kolmogorov (continuous line), Bohm (dashed line) and D(E) oc E 2 (dotted line). 
The separation between sources is d = 50 Mpc and the injection spectrum index is 7 3 = 2.7 
(see text). The AGASA-Akeno experimental data with the universal spectrum (dash-dotted 
line) are also reported. 



spectrum from E min ~ 1 GeV up to E^ x = 1 x 10 22 eV. For d = 50 Mpc L p = 1.5 x 10 49 erg/s 
and for d = 30 Mpc L p = 3.0 x 10 48 erg/s. To reduce these luminosities one can assume 
that the acceleration mechanism operates starting from some larger E m i n . Then the required 
luminosity is reduced by a factor E G ^f 2 \ which is 1.3 x 10^ 5 for E min = 1 x 10 8 GeV, and 
2.5 x 10~ 6 for E min = 1 x 10 9 GeV. Another possible assumption is the standard spectrum 
oc 1/E 2 at E < -E m i n as Berezinsky et al (2002b) have assumed. 

Figs 4 and 5 show that the dip seen in the universal spectrum as a signature of the in- 
teraction with CMB (Berezinsky et al 2002a) survives in the case of propagation in magnetic 
field with configuration (1 nG, 1 Mpc). As will be shown below the same is true for weaker 
and stronger magnetic fields. 

The case of a strong magnetic field (Bq, l c ) = (100 nG, 1 Mpc) is shown in Fig. 6. This 
is a very attractive case: the good agreement with the data is reached using the standard 
generation spectrum oc 1/E 2 and d = 100 Mpc. The required luminosity is reasonable, 
L p = 3 x 10 45 erg/s for E min ~ 1 GeV and E max = 1 x 10 22 eV. The diffusion coefficient used 
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Fig. 5. — The same as in Fig. 4, but for smaller separation between sources d = 30 Mpc. 
Notice the better conversion to the universal spectrum than in the case with d = 50 Mpc, 
shown in Fig. 4. 



in this case is D ps const at E < E c (the best fit in Fig. 6 is obtained for D(E) oc E om ). 
Unfortunately, the required magnetic field is much higher than that obtained in the MHD 
simulations by Dolag et al. (2004) and Sigl et al. (2004), though it does not contradict the 
existing observational upper limits. 

Let us now come over to the case of very weak magnetic field B ~ 0.1 nG, favored 
by MHD simulations by Dolag et al (2004). In this case E c « 1 x 10 17 (/ c /l Mpc) eV and 
kis(E) « 100E? 8 (1 Mpc// C ) Mpc. Therefore, for l c < 1 Mpc and E > 3 x 10 18 eV the protons 
propagate quasi-rectilinearly in the universe. In this case the distance between sources d is 
less than the propagation lengths ldm(E) and la,u(E), and the spectrum at least at energies 
(1 - 40) x 10 18 eV must be universal. 

A note of warning should be made about the validity of the Syrovatsky solution at 
E < 1 x 10 19 eV. This solution is expected to work not perfectly well at these energies, 
because it is valid only in the case when the energy losses b(E) and diffusion coefficient 
D(E) are time-independent 1 . For the above-mentioned energies this is not the case, because 



One should be careful with inserting ad hoc time-dependent quantities in the Syrovatsky solution (1), 
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Fig. 6. — Differential spectrum in the case of B = 100 nG, l c — 1 Mpc and for the diffusion 
regime with D(E) ~ const. The separation between sources is = 100 Mpc and the injection 
spectrum is the standard one (7 9 = 2). The AGASA-Akeno experimental data with the 
universal spectrum (dashed line) are also shown. The sharp cutoff in the energy spectrum 
at the highest energies is due to large distances to the nearby sources r ~ d ~ 100 Mpc. 

during the time of propagation the temperature of the CMB radiation changes appreciably, 
and hence the energy losses too 2 . The diffusion equation itself should be also modified as 
t — > to by the cosmological relations between time and distance. However, the approximate 
agreement, which we obtained (to be discussed somewhere else) between the Syrovatsky 
solution in quasi-rectilinear regime and the exact rectilinear propagation demonstrates the 
approximate validity of this solution at the discussed energies. 

Another argument in favor of the Syrovatsky solution as a reasonable approximation at 
the discussed energies E < 1 x 10 19 eV is the convergence to the universal spectrum (compare 
Fig. 4 and Fig. 5 for d = 50 Mpc and d = 30 Mpc, respectively). The universal spectrum is 



because in this case it ceases to be a solution of the corresponding diffusion equation. For example, it is 
forbidden to introduce the cosmological scaling factor a(t), because it results in time dependent energy losses 
b(E,t), or considering A in Eq. (1) as function of E,E g and t. 

2 In our previous paper, Aloisio and Berezinsky (2004), we deliberately limited ourselves to higher energies. 
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calculated in the case of time-dependent CMB temperature and for an expanding universe. 
The Syrovatsky solution converges to this spectrum with accuracy better then 15% when 
d — > 3 — 5 Mpc (to be discussed somewhere else). 

Following the papers by Berezinsky et al (2004) and Lemoine (2004), we shall now discuss 
shortly the transition from galactic to extragalactic cosmic rays. The remarkable feature of 
the diffusive spectra is the low-energy steepening at the fixed energy E s ~ 1 x 10 18 eV, which 
provide the transition from extragalactic to galactic CR. This energy coincides approximately 
with the position of the second knee E sk and gives a non-trivial explanation of its value as 

Like in the above-mentioned works we shall assume that at E > 1 x 10 17 eV the galactic 
spectrum is dominated by iron nuclei and calculate their flux by subtracting the calculated 
flux of extragalactic protons from all-particle Akeno spectrum. For these calculations we 
shall fix the spectrum with magnetic configuration (1 nG, 1 Mpc), the Bohm diffusion at 
E < E c and a separation between sources on the lattice d = 30 Mpc (see Fig. 5). The 
calculated spectrum of galactic iron is shown in Fig. 7 by the dashed curve. The fraction 
of iron-nuclei in the total flux is shown in Table 1 as a function of energy. This prediction 
should be taken with caution because of the model-dependent calculations (assumption of 
the Bohm diffusion) and uncertainties involved in the Syrovatsky solution. However, it is 
interesting to note that the iron-nuclei spectrum in Fig. 7 practically coincides with the 
spectrum calculated by Berezinsky et al. (2004) for the model with the generation spectrum 
steepening. The iron-nuclei spectra in both cases are well described by the Hall diffusion 
(Ptuskin et al. (1993)) in the galactic magnetic field at energies above the knee. 

We shall compare now our results with those obtained by Lemoine (2004), who also 
found the low-energy steepening of the spectrum due to diffusion. Lemoine has limited 
his calculations to the case Bq^/T c ~ 2 x 10 -10 GMpc 1 / 2 , while we demonstrated that this 
phenomenon is valid for much wider range of parameters, for example our configuration 
(100 nG, 1 Mpc) corresponds to the Lemoine parameter two order of magnitude larger. We 
considered here a more realistic basic scale l c ~ 1 Mpc and the various regimes of diffusion, 
while Lemoine limited himself to the D(E) oc E 2 regime only. We have also obtained the 
important result that the energy of the steepening is the same, E s ~ 1 x 10 18 eV, for all 
diffusion regimes and distances between the sources, and that universality is determined 
almost entirely by the proton energy losses. We discussed the diffusive anti-GZK effect, 
which we consider as the most interesting observation of this work 3 . 



3 Our work has been performed independently from the paper by Lemoine (2004) and much earlier. We 
discussed our results with Pasquale Blasi in September 2004. The delay with the publication was partly 
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\og[E (GeV)] 

Fig. 7. — The galactic iron-nuclei spectrum computed by subtracting the extralactic proton 
spectrum from the Akeno-AGASA data. The extragalactic proton spectrum is taken for the 
case B = 1 nG, l c — 1 Mpc, d = 30 Mpc, 7 9 = 2.7 with the Bohm diffusion at E < E c . 

5. Conclusions 

We have analyzed in this paper the anti-GZK effect in the diffusive propagation of ultra 
high energy protons. This effect consists in an increase of the maximum distance r max (E), 
from which ultra high energy protons can reach an observer, with an increasing of the energy 
E. This increase is terminated by a jump, which is located at energy Ej rs 2 x 10 18 eV. The 
position of the jump is determined exclusively by energy losses (transition from adiabatic 
to pair-production energy losses) and it is independent of the diffusion parameters. The 
position of the jump practically coincides with the position of the aforementioned transition 
and gives approximately the position of the second knee observed in the cosmic ray spectrum 
(see below). 

The observational consequences of the antiGZK effect is the low-energy "cutoff" of the 
diffuse spectrum, which is in fact a steepening in the spectrum, as the GZK cutoff is. The 
steepening energy E s coincides approximately with the position of the jump, E s ~ Ej, 
and it is also practically independent of the diffusion parameters, i.e. of the basic scale of 



connected with our attempts to overcome the problems of the Syrovatsky solution. 
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E (eV) 


10 17 


2 x 10 17 


5 x 10 17 


7 x 10 17 


10 18 


JFe/ (Jp + JFe) 


0.83 


0.66 


0.32 


0.17 


0.04 



Table 1: Fraction of iron-nuclei in the total flux as function of the energy. 

magnetic field coherence l c and of the magnetic field B on this scale. However, the shape 
of the steepening is determined by the diffusion regime: it is most steep in case D(E) oc E 2 
diffusion, most flat in case of the Kolmogorov diffusion, with the Bohm diffusion between 
them. 

In our calculations we have used the Syrovatsky solution to the diffusion equation, 
combined with the rectilinear propagation at the appropriate distances. The sources are 
located in the vertexes of a lattice with a spacing scale d (the source separation). We have 
used mostly d = 30 Mpc and d = 50 Mpc, which correspond to a space density of the sources 
3.7 x 10~ 5 Mpc -3 and 8.0 x 10~ 6 Mpc -3 , respectively. The observed small-angle clustering 
favors the density n s ~ (1 — 3) x 10~ 5 Mpc -3 . The diffusion coefficient D(E) is calculated 
for a random magnetic field with the basic scale l c and the coherent magnetic field on this 
scale Bq. Using this approach we have calculated the diffusive spectra for various magnetic 
configurations (B , l c ) and source separations d. 

Physically the most reasonable case corresponds to a magnetic field configuration (1 nG, 
1 Mpc) with a source separation d = 30 Mpc and d = 50 Mpc. The calculated spectra are 
shown in Figs 4 and 5 in comparison with Akeno-AGASA data. For a power-law generation 
spectrum with 7 S = 2.7 the agreement is good, but needs too high luminosity of the sources 
L p , if the power-law spectrum starts with low energy E min ~ 1 GeV. This problem can be 
amiliorated assuming higher values of E min . 

The calculated diffusive spectra in the energy interval (1 — 80) x 10 18 eV agree perfectly 
well with the universal spectrum and experimental data, showing the presence of the dip 
caused by e + e~ production. 

An interesting case is given by the diffusion in strong magnetic field with basic config- 
uration (100 nG, 1 Mpc) and source separation d=100 Mpc. In this case (Fig 6) the best 
fit of the spectrum is obtained for the standard acceleration spectrum Q(E) oc 1/E 2 and 
E min ~ 1 GeV. The required luminosity is L p = 3 x 10 45 erg/s. Up to energy E ~ 1 x 10 20 eV 
the predicted spectrum agrees with data of both detectors, AGASA and HiRes. The sharp 
cutoff at E ~ 1 x 10 20 eV is produced due to large distances r ~ d to the nearby sources. 
For the explanation of the AGASA excess at E > 1 x 10 20 eV a new component of ultra high 
energy cosmic rays (e.g. from superheavy dark matter, see Aloisio et al 2004) is needed. 
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At energies E < E s , where Uis becomes much smaller than d, the diffusive spectrum 
exhibits a steepening in contrast to the universal spectrum (see Figs 4-6). 

The steepening of the spectrum at E s ~ 1 x 10 18 eV provides a natural transition from 
galactic to extragalactic cosmic rays. This energy coincides with the second knee observed 
in cosmic rays spectra by most of the detectors. While the energy of the transition E s (and 
thus position of the second knee) is predicted in a model independent way, the shape of 
the proton spectrum below 1 x 10 18 eV and the fraction of galactic iron nuclei are model 
dependent: they differ for various diffusion regimes. 
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